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ABSTRACT 

We investigate how the probability of the formation of giant arcs in galaxy clusters is expected to change in 
cosmological models dominated by dark energy with an equation of state p — wpc 2 compared to cosmological- 
constant or open models. To do so, we use a simple analytic model for arc cross sections based on the Navarro- 
Frenk- White density profile which we demonstrate reproduces essential features of numerically determined arc 
cross sections. Since analytic lens models are known to be inadequate for accurate absolute quantifications of arc 
probabilities, we use them only for studying changes relative to cosmological-constant models. Our main results 
are (1) the order of magnitude difference between the arc probabilities in low density, spatially flat and open CDM 
models found numerically is reproduced by our analytic model, and (2) dark-energy cosmologies with w > — 1 
increase the arc optical depth by at most a factor of two and are thus unlikely to reconcile arc statistics with 
spatially flat cosmological models with low matter density. 



1. Introduction 
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Q_The statistics of giant luminous arcs in the cores of galaxy clus- 
1. ters has long been recognised as a po t entially powerful cosmo - 
^ logical probe (e.g. Wu & Maolll996t iBartelma nn et alJ ll998). 

Arcs are formed by gravitational lensing from sources which 
^ happen to lie close to the caustic curves of a cluster lens, where 
• • magnification and distortion are particularly strong. 
. ' The cosmological power of arc statistics derives from at least 
k> two, maybe three principal reasons. First, for clusters to be ef- 
5_j ficient lenses, they have to be located approximately half-way 
C$ between the sources, typically around redshift unity, and the 
observer. Depending mostly on the mean cosmic matter den- 
sity, parameterised by £2o, clusters form earlier or later in cos- 
mic history if the matter density is low or high, respectively. In 
high-density model universes, the cluster population at the red- 
shifts mostly relevant for lensing, z ~ 0.3 — 0.4 is substantially 
less rich than in low-density universes, reduci ng the number 
of available efficient lenses dramatically (e.g. iRichstone et alJ 
[T992tlBartelmann et alJl993linc7v & Colell 993111 9941) . 

The second principal reason is that strong lensing is a highly 
nonlinear phenomenon in the sense that it requires the lensing 
mass distribution to be supercritical for strong lensing, which 
means that a suitable combination of surface mass density and 
gravitational tidal field needs to be large enough, and that, once 
a lens is supercritical, even small changes in both can change 
significantly the length of the caustic curves, and thus the lens' 
ability for strong lensing. 

Different cosmological models predict the mass distribution 
in clusters to be more or less concentrated. Numerical sim- 
ulations consistently show that, the earlier a dark-matter halo 
forms, the more concentrated it is because it appears to keep a 
record of the mean cosmic density at the time when it formed 
JNavarro et all 119961 119971) . Structure forms later in spatially 
flat than in open, low density cosmological models, thus halos 
in models with cosmological constant are generally less concen- 



trated than halos in open models. 

A possible third reason is that the gravitational tidal field 
at the location of the lens plays a very im portant role 
JBartelmann et all I1995L iMeneghetti et all 120011) . It is strong 
if lenses are highly asymmetric, as clusters frequently are, and 
if the surrounding matter distribution is highly irregular. It is 
possible that, if cosmic structure forms later, cluster mass distri- 
butions are less relaxed and thus more asymmetric, and that also 
the degree of irregularity in their neighbourhood is different than 
if structures formed earlier. On the other hand, clusters forming 
earlier are built from subhalos which tend to be more concen- 
trated and thus more strongly gravitationally bound, hence sub- 
structures could then persist within clusters for a longer time and 
contribute to the asymmetry. 

Unfortunately, the combination of these effects renders an- 
alytic models for arc statistics entirely inade quate for accurate 
quantitative predictions of arc probabilities (Meneghetti et al., 
2003). The effects of cosmology on cluster compactness and 
asymmetry, and on the tidal field of the matter surrounding the 
clusters, cannot be captured by reasonably simple analytic lens 
models. Numerical simulations of arc statistics, using clusters 
formed in sufficiently large A^-body simulations as lenses, led to 
the surprising result that the expected number of giant luminous 
arcs on the sky differs by orders of magnitude between different 
cosmological models. While a model with critical matter density 
and no cosmological constant fell below the observed number of 
arcs, extrapolated to the full sky, by two orders of magnitude, a 
ACDM model with £2o = 0.3 failed by one order of magnitude, 
and only a low-density open model with £l p = 0.3 produced ap - 
proximately the right number of arcs IBartelma nn et all ll998). 

The statistics of quasars multiply imaged by galaxies has of- 
ten been used for constraining cosmological parameters. The 
basic argument is that the number of lenses and their redshifts 
should increase as Ha increase s, which typically yields up- 
per limits on Ha 5, 0-6 — .7 (Kochanek, 119961 iFalco et all 
1998; Ouast & Helbig lll999). although discrepant results have 



also b een found (iChiba & Yoshiil H999t iHelbid U999t iKeetonl 
2002). We emphasise that the sensitivity of cluster lensing to 
is of a different nature. Since clusters form much later in cosmic 
history than galaxies, the volume effect is negligible, but Q.\ 
changes the dynamics of cluster formation and thus their core 
structure, to which strong lensing is highly sensitive. 

These numerical res ults of [B artelmann e t all (119981) w ere 
tested by Coorav ( 1999) and Kaufmann & Straumann ( 2000) us- 
ing analytic models based on singular isothermal spheres. They 
could confirm the sensitivity of arc statistics to £2o, but found 
only a very weak dependence on Ha, in contrast to the numer- 
ical results. The isothermal sphere has two disadvantages with 
respect to arc statistics. First, arc cross sections are very sensi- 
tive to asymmetries in the cluster mass distribution, thus axially 
symmetric models lack a property which is crucially important 
for arc statistics. Second, numerical simulations show that the 
central density concentration of clusters depends on cosmology, 
and this potentially important feature is not reproduced by the 
scale-free isothermal models either. We shall construct in this 
paper an analytic model which qualitatively reproduces the ear- 
lier numerical results. 

iBartelmann et alJ {1298) used two completely different types 
of A^-body codes for simulating galaxy clusters. Numerous 
subsequent tests of the results showed that the arc numbers 
derived could be off by factors of perhaps 1 .5 to 2, but that 
there was n o way how order -of-ma gnitude diffe rences could 
be bridged JMeneghetti et all 1200(1 iFlores et all 120001) . The 
problem became substantially more acute when measurements 
of the cosmic microwave background (CMB), combined with 
observations of supernovae of type la and large-scale galaxy 
surveys, left very little room for model universes which are 
not spatially fiat and hav e density parameters much d iffer- 
ent from Q.t) ^ 0.3 (e.g. iRiess et alJ Il998t iPerlmutter et al] 
[19991 iLee et alJpOOlttAbroe et al 1120021: lEfstathiou et all2 002: 
Nett erfield et al.l2002t1Wang et alJ2002TlHalverson et all2002l) 
Obviously, there is an interesting discrepancy between the statis- 
tics of arcs seen on the sky, and the probability for arcs pro- 
duced in cosmological models which are convincingly required 
by various other observations. While observations consis- 
tentl y indicate a high probability for arc formation in clusters 
*Le Fevre et all H994tlGioia & Lupmno, 1994; Luprjino et all 



1999; IZaritskv & Gonza lez. 2003; Gladd ers et all |2003), we 
should point out that the discrepancy between theory and obser- 
vations is so far only based on a single set of simulated clusters. 

A spatially fiat universe with low matter density and a cos- 
mological constant is extremely difficult to justify theoretically. 
The vacuum energy density provided by the cosmological con- 
stant is tens of orders of magnitudes below any natural scale 
which is conceivable in particle physics (see |Carroll 2002 f° r 
a review on the cosmological constant problem). This diffi- 
culty motivated the introduction of a more general concept for 
a vacuum energy cosmological component, now widely known 
as dark energy. The theoretical and observational aspects of the 
dark energy are one of the most im portant issues in modern cos- 
mology (see Peebles & Ratra 2002 for a review). In general, 
the most important difference of a dark energy component com- 
pared to a cosmological constant is that its equation of state, 
w, can be different from —1, generally implying a time vari- 
ation. It should be noted that recent analyses of CMB data 
seem to favour a value of w very close to —1, albeit these re- 
sults suffer to some degree from parameter degen eracies an d are 
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Snergel et all 120031). 


Recently, Bartelmann et al. 


( 2002) argued that dark matter ha- 



ergy density today. The main reason is that halos form earlier in 
dark energy models, allowing them to be more compact. Sev- 
eral different recipes for describing halo concentrations found 
in numerical simulations as a function of their formation time 
lead to consistent results. The halo concentration increases no- 
ticeably in the interval — 1 < w < —0.6. For higher values of 
w, which ar e too high for the cosmic acceleration to agree with 
recent data (Riess et al., 1998; Perlmutter et al., 1999), there is a 
strongly opposing effect related to the amplitude of fluctuations 
in the CMB: The high level of the Integrated Sachs-Wolfe (ISW) 
effect on the large scale CMB anisotropies leads to a sharp de- 
crease in the norm alisation of the dark matter power spectrum 
( Bartelma nn et all 120021) . Within a cosmologically interesting 
range for the equations of state of the dark energy, the balance 
between the ISW effect and the earlier formation of halos is 
quite delicate, but halos can typically be expected to be more 
concentrated. 

The ability of a galaxy cluster to produce giant arcs depends 
sensitively on the concentration of its mass profile because of 
the nonlinearity of the strong lensing effect. We therefore wish 
to investigate how the probability for arc formation changes in 
dark energy models, compared to cosmological-constant or open 
models. For doing so, we use a simple, analytic description 
for the arc cross section of a cluster of given mass, which we 
demonstrate to possess the relevant features of the fully numer- 
ical results. Although it has been shown that analytic models 
are inadequate for quantitatively reliable arc statistics, we are 
here interested only in the relative change of the arc-formation 
probability caused by changes in the cosmological model. 

Section 2 of the paper introduces the cosmological back- 
ground model. Section 3 describes our simple analytic model 
for the arc cross section of a galaxy cluster. We then use this 
model in Sect. 4 for calculating arc probabilities, and summarise 
our conclusions in Sect. 5 

2. Cosmological Model 

Dark energy is characterised by a negative pressure, p = wpc 2 , 
where pc 2 is the mean energy density of the universe and the 
equation of state w assumes negative values in order to produce 
cosmic acceleration according to the data from type la super- 
novae (Riess etal., 1998; Perlmutter et al., 1999). Theoretical 
models of dark energy, such as Quintessence scalar fields, gen- 
erally predict a time variation of the equation of state, as well as 
the presence of dark energy fluctuatio ns on supe r-horizon cos- 
mological scales (see e.g. Peebles & Ratra 2002 and references 
therein). In this work, we concentrate on the very basic aspect 
of dark energy, by neglecting the spatial inhomogeneities and 
assuming w to be a constant. In this case, the adiabatic equation 
requires the equivalent matter density Pq of the dark energy to 
change with the cosmological scale factor a as 



PQ = PQ.oa 
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starting from the density Pq o today. Obviously, cosmological 
constant models are retained setting w = — 1. Replacing the 
conventional cosmological-constant term by a dark-energy term, 
Friedmann's equation reads 



H 2 {a)=Hl 
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los in simple dark energy models should be more concentrated 
than in cosmological constant models with the same dark en- 



assuming with support from recent measurements of 
anisotropies in the cosmic microwave background that the 
curvature term is negligible. Here, H(a) is the Hubble parame- 
ter as a function of a, Ho is the Hubble constant, and £2o is the 
density parameter for non-relativistic matter. 



The main consequences for the structure formation process of 
this modified term in Fried mann's eq uation have been detailed 
in an earlier paper (Bartelmann et al., 2002), so we summarise 
them only briefly here. Starting from a low-density cosmologi- 
cal constant model, keeping fixed and increasing w, the cos- 
mic volume per unit redshift shrinks, and the linear growth factor 
for cosmic structures starts rising earlier, hence structures start 
forming earlier if the cosmological constant is replaced by a dark 
energy component. Additional effects are that the parameters 
characterising halo formation, i.e. the linea r overdensity S r and 
the v i rial overdensity A v , are changed (cf. IWang & Steinhardt 
ll998HMainini et all2003h . 

As we anticipated in the introduction, a further important con- 
sequence of dark energy models is that the gravitational poten- 
tial of the density fluctuations changes more rapidly with time, 
leading to an increased integrated Sachs-Wolfe effect on the 
CMB fluctuations at large angular scales. Given the observed 
level of CMB power on large angular scales 1 , a decreasing frac- 
tion of the observed anisotropies can thus be attributed to the 
primordial CMB fluctuations, hence the normalisation of the 
power spectrum has to be reduced as w is increased in order to 
be compliant with the COBE-DMR data. The main result is that 
the power-spectrum normalisation Os is decreasing gently as w 
is increased from — 1 to ~ —0.6, and then turns to drop more 
rapidly as the ISW effect intensifies. Here and below, we adopt 
the CDM power spectrum with the Harrison-Zer dovich power- 
law in dex n — 1 and the transfer function given by Bardeen et al. 
( 1986). Throughout, we use Q.q = 0.3 and £2q = 0.7. It should 
be noted here that our results will be sensitive to the exact value 
of n as well as other cosmological parameters. 

3. Arc Cross Sections 

3.1. Halo Model 

We will assume in the following that the average radial density 
profile of galaxy clusters can be described by th e profile found in 
numerical simulations by Navarro et alJ i 19961 hereafter NFW), 



(3) 



r/r s (l+r/r s ) 2 ' 



where p s is a density scale, and r s a scale radius. The ratio be- 
tween r s and the radius r2oo enclosing a mean halo density of 
200 times the critical density is called concentration, c — r & /ryxi- 
The two parameters, p s and r s , are not independent. Numerical 
simulations showed that the concentration parameter c depends 
on the halo mass, which is thus the only free parameter. 

Numerically simulated halos tend to be the more concentrated 
the earlier they form. Their central density apparently reflects 
the mean cosmic density at the time of their formation. Since 
halos of higher mass are formed later in hierarchical models than 
halos of lower mass, the concentration is decreasing with halo 
mass. Based on these findings, several algorithms were designed 
for computing halo concentra tions from h alo masses. 

The algorithm by Navarro et al. ( 1997) first assigns to a halo 
of mass M a collapse redshift z co n defined as the redshift at 
which half of the final halo mass is contained in progenitors 
more massive than a fraction /nfw of the final mass. Then, the 
density scale of the halo is assumed to be some factor C times the 
mean cosmic density at the collapse redshift. They recommend 
setting /nfw = 0.01 and C — 3 x 10 3 because their numerically 

'The only existing data are from the Differential Microwave Radiometer 
(DMR) on board the COsmic Background Explorer satellite (COBE, see Smoot 
1999 and references therein); more data will be available in the near future from 
the Microwave Anisotropy Probe (MAR map . gsf c . nasa . gov) and the Planck 
(astro . estec . esa. nl/SA-general/Pro j ects/Planck) satellites. 



determined halo concentrations were well fit assuming these val- 
u es. 

iBullock et al.l (120011) suggested a somewhat simpler algo- 
rithm. Haloes are assigned a collapse redshift defined such that 
the non-linear mass scale at that redshift is a fraction /b of the 
final halo mass. The halo concentration is then assumed to be a 
factor K times the ratio of the scale factors at the redshift when 
the halo is identified and at the collapse redshift. Comparing 
with numerical simulations, they found /b = 0.01 and K = 4. 
This algorithm reflects the change of halo concentration s with 
redshift more accurately than the approach by iNavarro et alJ 
ill 9971) predicts. 

A third algorithm was suggested bv lEke et alJ J2001I) . They 
assigned the collapse redshift to a halo of mass M by requir- 
ing that the suitably defined amplitude of the linearly evolving 
power spectrum at the mass scale M equals a constant C^ s . 
Numerical results are well represented setting Cens = 28. 
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FIG. 1. — Halo concentrations according to Eke et al. (2001) are shown 
as functions of halo mass for four spatially flat cosmological mod- 
els with different choices for the dark energy parameter w and for 
the OCDM model, as indicated. As w increases, halos become 
more concentrated until w ~ —0.6. If w increases further, halo con- 
centrations drop because then the amplitude rjg of the power spec- 
trum has to decrease rapidly in order to remain consistent with the 
COBE-DMR data, as the integrated Sachs-Wolfe effect becomes larger 
(cf. lBartelmann et all2002T) . 

We will adopt the latter method for this paper because the con- 
centrations computed from the algorithm by Navarro et al. drop 
too slowly with redshift compared to numerical simulations, 
and the algorithm by Bullock et al. has problems for high halo 
masses because of the requirement that a fixed fraction of the fi- 
nal halo mass should equal the nonlinear mass, which may never 
be reached if the halo mass is high. Concentrations as a func- 
tion of halo mass for four different choices of the dark energy 
parameter w are shown in Fig.^ 



3.2. Strong Lensing by NFW Haloes 

The arc cross section of a galaxy cluster is the area in the source 
plane where a source has to lie for being imaged as an arc with 
specified properties, e.g. exceeding a threshold length-to-width 
ratio. A typical cluster lens has two critical curves, defined as 
curves in the lens plane along which the Jacobian matrix of the 
lens mapping is singular, and the image magnification is for- 
mally infinite (cf. Fig.|2j. The critical curves are the images of 
the caustic curves, thus sources close to a caustic are imaged as 



highly magnified and distorted images. The two critical curves 
and their corresponding caustic are called tangential and radial, 
because of the dominant orientation of the image distortion rela- 
tive to the centre of the lens. Seeking to quantify large-arc cross 
sections, we are thus looking for an appropriately defined area 
covering the tangential caustic of a cluster. 

It is well known that arc cross sections depend strongly on the 
exact shape of the cluster mass distribution, and on the gravita- 
tional tidal field exerted by de nsity fluctuations in its neighbour- 
hood (Bartelmann et al., 1995). Thus, the only reliable method 
for exactly quantifying arc cross sections has to use numerically 
simulated cluster models without referring to any symmetry as- 
sumptions. For a study like ours, however, we only need to de- 
scribe how arc cross sections are expected to change relative to a 
fiducial model when certain cosmological parameters are modi- 
fied. 
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FIG. 2. — Tangential and radial caustics for NFW lens models with el- 
liptically distorted lensing potential. A halo of mass 10 15 ^" 1 M Q at 
redshift z = 0.3 is assumed here, and the sources are all placed at z s = 1 . 
The cosmological model is spatially flat with Q.q = 0.3 and normalised 
to the COBE-DMR data. The upper panels show the caustics for a 
dark energy parameter of w = — 1 and for two ellipticities, e = 0.1 and 
e = 0.4. The lower two panels show the caustics for w = —0.6. The 
figure illustrates the high sensitivity of strong lensing to halo ellipticity 
and concentration. 

Gravitational lensing b y an NFW halo can be described by 
its le nsing potential (e.g. Bartelmanr] 1 1 99(1 iMeneghetti et alJ 
l200l . 



\\f(x) = 4k s 
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which is related to the lensing convergence through the Poisson 
equation 



K(x) = 

The factor K s in Eq. 10} is defined by 

Ps^s 



(5) 



(6) 



where E cr is the critical surface mass density for lensing. The di- 
mensionless radius x = r/r s can conveniently be replaced by the 
angular radius = r/D^ = xr^/D^, where D& is the angular di- 
ameter distance from the observer to the lens. The gravitational 
tidal field, or shear, of the lens is the two-component quantity 



Yi = ^(¥n -¥22) , Y2 = ¥i2, 



(7) 



where the subscripts abbreviate partial derivatives with respect 
to the angular coordinates (01,82) on the sky. The deflection 
angle is the gradient of the lensing potential, 



a(8) = V\|/(0) . 



(8) 



The Jacobian matrix of the lens mapping has the components 



% = 8y-¥i; 



(9) 



Its eigenvalues are A,± = (1 — k) ±y, where Y= (Yi +^) 1,/2 i s 
the amplitude of the shear. The tangential critical curve is deter- 
mined by the condition X_ = (1 — k) — y = 0. 

We now distort the axially symmetric NFW lens such that the 
iso-potential lines become ellipses, 
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As noted by iKassiola & Kovne'rl (119931) . an elliptical poten- 
tial can lead to dumbbell-shaped mass distributions with lo- 
cally negative mass density. As real clusters are irregular, 
dumbbell-shaped mass distributions are acceptable. For the 
NFW profile with elliptical isopotential contours, the mass den- 
sity does indeed become mildly negative, but only well outside 
the core where strong lensing occurs. For e = 0.4, the mini- 
mum K is ~ —0.01 t imes the convergence in the core (see also 
Idolse & K neib 2002). We thus use the elliptical lensing poten- 
tial (II 01 for computational simplicity. 

The Jacobian matrix and its eigenvalues can be computed 
from dlOt using the relations introduced before. Generally, the 
zeroes of the tangential eigenvalue A,_ have to be determined 
numerically. On the coordinate axes, they are given by 



f (l-e)oc (02=0) 
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The corresponding caustic points, i.e. the cusps of the diamond- 
shaped caustic on coordinate axes, can then be found using the 
lens equation, 

P = 0-a(0). (12) 

We thus know the four intersection points of the tangential caus- 
tic curve with the coordinate axes. We defined them to lie at 
(01,02) = (0,±fl) and (0i,02) = (±b,0). Since the major axis 
of the iso-potential ellipses points along the 02 axis, the cusps 
on the 02 axis are further away from the lens centre than the 
cusps on the 0i axis, hence a> b. Figure|2]illustrates the caustic 
curves for elliptical NFW lens models with two different ellip- 
ticities in two spatially flat cosmological models with different 
values for w. 

A simple assumption for the large-arc cross section a of the 
elliptically distorted NFW lens holds that it is proportional to 
the area enclosed by the critical curve, hence a °c ab, with a pro- 
portionality constant depending on the exact shape of the caus- 
tic curve, and thus on the ellipticity of the lens model. Since 
we do not require an absolute calibration of the arc cross sec- 
tions, we adopt o = ab. Earlier work ( Menegh etti et all [2000: 



iFlores etaQ 120001) has shown that individual galaxies have a 
negligible effect on arc cross sections, which further supports 
the assumption that they are determined by the overall extent of 
the caustic curves. 

We will later have to integrate over the cluster population, 
conveniently parameterised by the virial mass and described by 
the mass function. Thus, we have to verify whether our approxi- 
mate description of the arc cross section scales with cluster mass 
in the same way as numerically determined arc cross sections of 
the same cluster models, given identical ellipticity parameters 
e. We therefore set up deflection-angle maps starting from the 
elliptically distorted lensing potential \(/($) and used them for 
imaging randomly distributed, intrinsically elliptical sources by 
tracing rays passing the image plane on a rectangular grid. The 
sources are placed on adaptively refined grids whose resolution 
is progressively increased in the vicinity of caustic curves. The 
images are automatically classified according to their length, 
width, length-to-width ratio and several other parameters, and 
the cross section is determined by counting the number of im- 
ages exceeding a threshold length-to-widt h ratio. More detail 
on this method can be found in (Bartelma nn & WeissL H9941 
iMeneehetti et all l2000). Results are shown in Fig. [3] 
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FIG. 3. — Comparison between the numerically determined arc cross 
sections of elliptically distorted NFW lenses as a function of their virial 
mass M (solid curve), and the simple estimate for the arc cross sections 
which is proportional to the area covered by the tangential caustic curve 
(dotted curve). Both curves are arbitrarily normalised to unity at the 
high-mass end. The lenses are placed at redshift z = 0.3, the sources 
at redshift z s = 1.0, the ellipticity of the NFW potential is e = 0.4. 
The curves were computed for a dark energy cosmological model with 
w = —0.6. Evidently, the simple estimate for the cross sections cor- 
rectly reproduces the scaling of the cross sections with cluster mass. 
The numerically determined curve is not smooth because the cross sec- 
tions are computed from finite numbers of simulated arcs and thus sub- 
ject to random fluctuations. 

The essential feature of the numerically determined arc cross 
sections is that they increase approximately quadratically with 
the lens mass. Since the Einstein radius of an extended lens 
typically scales linearly with the lens mass, this indicates that 
the arc cross section scales approximately like the square of the 
Einstein radius. Our simple estimate for the arc cross section is 
defined to reproduce this property. 

We also check whether numerically determined cross sec- 
tions scale with the dark energy parameter w in a similar way 
as our simple cross-section estimate does. Figure 0] shows the 
example of an elliptically distorted NFW lens of mass M = 
7.5 x 10 14 /i _l M©. The curves plotted there show qualitatively 



the same behaviour. As w increases above —1, the arc cross sec- 
tion increases by a factor of ~ 1 .9 until w rj —0.6, and then drops 
as the normalisation constraint requires to significantly reduce 
Og. Both curves are arbitrarily normalised to unity at their peaks. 
The differences between the numerically determined cross sec- 
tions and our simple estimate are unlikely to be significant be- 
cause the numerical method uses the images of a finite num- 
ber of discrete sources for determining the cross sections, which 
causes noise in the results. In any case, the two essential prop- 
erties of the change in the cross sections with w are well re- 
produced, namely the amplitude of the relative increase and the 
location of the peak. 
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FIG. 4. — Comparison between cross sections similar to Fig. [5] but for a 
fixed halo mass of M= 7.5 x 10 14 /i _1 Mq and varying dark energy pa- 
rameter w. The solid curve showing the numerically determined cross 
sections closely follows the dotted curve, which represents the simple 
cross-section estimate introduced here. The curves are arbitrarily nor- 
malised to unity at their starting point, i.e. at w = — 1. Increasing w 
from —1 to —0.6 increases the cross sections by a factor of ~ 1.9. As 
in Fig. [3] lens and source redshifts are set to 0.3 and 1, respectively, and 
the ellipticity of the NFW lensing potential is set to e = 0.2. 

Finally, we set the cross section to zero if the major axis of the 
critical curve falls below some threshold 8 m ; n . The idea behind 
this is that if the critical curves become too small, the images 
near the critical curves can hardly be called giant arcs. Suppose 
typical sources have diameters on the order of an arc second, 
and the lens should be able to produce arcs with a length-to- 
width ratio around ten. Then, ignoring the source magnifica- 
tion in the radial direction, the tangential critical curve needs to 
have a radius of approximately 10" for this to happen. We thus 
set m ; n = 10" unless stated otherwise, and show the effect of 
changing m i n to 5" below. 

We conclude from this section that the scaling of our sim- 
ple estimate for arc cross sections with lens mass and with the 
dark energy parameter w well reproduces what is expected from 
numerical treatments of the same lens models, i.e. NFW lenses 
with elliptically distorted lensing potential. We emphasise again 
that the absolute value of the cross sections are unimportant for 
our present purposes, as we are aiming at studying the change 
in the arc-formation probability in various cosmological models 
relative to the ACDM model. 

4. Arc Probabilities 

We can now proceed to compute the probability for arc for- 
mation by a population of clusters. Given a mass function 



dn(M,z)/dM, we can write the so-called optical depth x(z s ) as 



dz(l 



dV 
dF 



A/,„ 



dM^a(M,z), 



(13) 



where z s is the source redshift, V is the cosmic volume, and 
the factor (1 + z) 3 accounts for the fact that the mass function 
is defined per comoving volume. The lower mass limit M m j n 
is determined by the mass required to produce critical curves 
whose major axis exceeds the threshold 9 m ; n introduced in the 
preceding section, thus it depends on the lens redshift z. 

For the mass function dn / d M, we choose the modificat ion by 
Shet h & Tormer] i 1999ft of the iPress & Schechterl i 19741) mass 
function. It well reproduces the halo mass functions found in 
numerical simulations. We take into account that our defini- 
tion of mass differs slightly from Sheth & Tormen's in that we 
use the mass enclosed by a sphere in which the mean density is 
200 times the critical rather than the mean density. The mass 
function depends on cosmology through the normalisation of 
the power spectrum and the linear overdensity parameter derived 
from the spherical collapse model. 

We note an important difference to strong lensing by galax- 
ies. While the population of galaxy lenses is well de- 
scribed as isothermal sphe res with number counts de rived 
from observations (e.g. iLehar et alJ feOOft iKeetonl 12002). the 
lack of wide-separation lenses argues against isothermal den- 
sity profiles in cluster-s cale lenses t Flore s & Prim ack, 1994; 
Porciani & Madau, 2000). Baryonic physics changes the central 
density profiles of galaxy-scale halos, b ut is inefficient on clus- 
ter scales (e.g. Kochanek & Whita 120011) . Cluster mass func- 
tions derived from X-ray observations are found to agree well 
with theoretical predictions based on Press-Schechter type mod- 
els (e.g. Reiprich & Bohringer 2002). Thus modelling the clus- 
ter population with NFW density profiles and the Sheth-Tormen 
mass functions is well justified. 

We show in Fig.[5]the redshift integrand of Eq. dl3> . i.e. the 
cosmic volume times the integral over mass of the mass function 
times the arc cross section. The figure thus illustrates the total 
arc cross section contributed by the cluster population at redshift 
Z. For simplicity, sources are assumed to be at a single redshift of 
z s = 1 here. The four curves in Fig. [5] are for the ACDM model 
with CIq — 0.3, the open CDM model with the same D.q but with 
£Ia = 0, and for dark energy models with the same £2o but CIq — 
0.7, and w = —0.6 or w = —0.8. Again, the curves are arbitrarily 
normalised such that the optical depth x of the OCDM model is 
unity. The value of G% = 0.88 of the ACDM model was also set 
for the OCDM model for easier comparison. The COBE-DMR 
data would require as ~ 0.4 for the OCDM model, which is way 
below the value required for reproducing the observed number 
density of massive clusters. We thus have to choose Os for the 
OCDM model in conflict with the COBE-DMR data in order to 
produce comparable results on arc statistics. 

Figure[5]shows two important results. First, the simple model 
for arc cross sections introduced here is capable of reproducing 
the order-of-magnitude difference in the total arc cross section 
between the ACDM and the OCDM models that had been found 
earlier in numerical simulations, and could not be reproduced 
by analytic models based on singular isothermal cluster mass 
distributions. Second, although the dark energy models have a 
somewhat higher total arc cross section than the ACDM model, 
they are still by a factor of ~ 6 below the arc cross section for the 
OCDM model. According to our analytic estimates, dark energy 
models are thus unable to reconcile spatially flat cosmological 
models with low matter density with the high abundance of large 
arcs, which seems to be similarly impossible with cosmological 
constant models given the earlier numerical results. 




0.1 



0.01 



FIG. 5. — The redshift integrand of Eq. Jl 3 1 is plotted as a function of z 
for four different cosmological models, a ACDM model with iio = 0.3, 
and open CDM model with £Iq = 0.3 and no cosmological constant, and 
two spatially-flat dark energy models with Q.q = 0.3 and w = —0.8 and 
w — —0.6. The spatially-flat models are normalised to the COBE-DMR 
data, and the open CDM model has the same Og as the cosmological- 
constant model for easier comparison. The curves show that our simple 
analytic model succeeds in reproducing the order-of-magnitude differ- 
ence between the open and the cosmological-constant model found in 
numerical simulations, and that spatially flat dark energy models cannot 
bridge the gap between these two models. The ellipticity of the NFW 
lensing potential was set to e = 0.3 here. 



Figure |S] shows the optical depth x as defined in Eq. ( 1131 for 
four different choices of the free parameters we have introduced, 
namely the ellipticity e of the lenses and the cutoff radius 8 m i n . 
All curves have in common that the arc optical depth increases 
noticeably as w increases from — 1 to ~ —0.6, and drops rapidly 
as w is increased further. The curves are intended to show the 
relative increase in x compared to the ACDM model and are thus 
normalised to unity at w = — 1 . The peak amplitudes range from 
1.5 to 2. The largest increase is achieved for 9 m i n = 10" and 
e = 0.3, the smallest for 8 m ; n = 5" and e = 0.4. 

Several effects act jointly here. First, we saw in Fig. 0] that 
the arc cross section of (one example for) an individual halo in- 
creases by a factor of < 2 as w is increased from —1 to —0.6. 
An additional effect is that, as the halo concentration increases, 
halos of lower mass become capable of strong lensing. Since 
the mass function of galaxy clusters is steep, a small extension 
of the mass range towards lower masses can markedly increase 
the number of clusters available for strong lensing, but the re- 
quirement that arcs should be large imposes a lower limit on the 
cluster masses. We see the combined effect in Figs.|5]and|6] If 
we set the cutoff radius to 8 m i n = 10", we select for higher-mass 
clusters in the first place, whose mass function is steeper than for 
lower-mass clusters. Thus, the effect of lowering the lower mass 
limit by increasing the halo concentrations is more pronounced 
if the cutoff radius is chosen higher. 

Increasing the ellipticity of the lensing potential increases its 
gravitational tidal field, or shear. The increase in halo concentra- 
tion caused by the earlier halo formation in dark energy models 
with w > — 1 is then relatively less important for arc formation, 
which explains why the arc optical depth caused by cluster pop- 
ulations with lower e changes more strongly with w than for 
clusters with higher e. Fitting the elliptical NFW model to nu- 
merically simulated clusters yields values for e closer to 0.3. 
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FIG. 6. — Optical depth x(z s = 1) defined by Eq. (HJ for spatially-flat 
dark energy models as functions of w. All curves are arbitrarily nor- 
malised to their values for w = — 1, i.e. they show the change in optical 
depth with w relative to the cosmological-constant model. All curves 
were obtained for £Iq = 0.3 and COBE-normalised CDM power spec- 
tra. Results are shown for two different choices each for the ellipticity 
e and the cutoff angle 6 m ; n , as indicated. The curves show that clusters 
in dark energy models can be ~ 50% up to ~ 100% more efficient in 
forming arcs than in cosmological-constant models, depending in detail 
on the exact choices for the ellipticity and the cutoff angle. 



5. Summary and Conclusions 

We investigated how the probability for the formation of large 
gravitational arcs in galaxy clusters is expected to change as the 
underlying cosmological model is modified. The main reason 
for this investigation was our earlier finding that halo concentra- 
tions in a simple class of spatially flat, dark-energy dominated 
cosmological models are expe cted to depend sensitivel y on the 
equation-of-state parameter w (Bartelmann et al., 2002). Strong 
lensing in general, and the formation of large arcs in particular, is 
a highly nonlinear effect which depends sensitively on the matter 
concentration in the lens cores. It could therefore reasonably be 
expected that arc probabilities would change significantly with 
w. 

This is o f cosmological re l evance since earlier numerical 
simulations (Bartelmann et al., 1998) showed that clusters in 
cosmological-constant models fall short by about an order of 
magnitude of reproducing the abundance of observed large grav- 
itational arcs. While virtually all recent cosmological exper- 
iments favour a spatially flat, low-density universe, arc statis- 
tics apparently strongly prefers a low-density, open model over 
a low-density model with cosmological constant. This nu- 
merical result was questioned on the basis of simple analytic 
models which failed to reproduce the strong dependence of 
arc statistics on Qa found in the simulations (iCooravl [1999; 
Kaufmann & Straumannf 120001) . Our first goal was thus to in- 
vestigate whether the numerical results could be supported by a 
more detailed analytic model. 

Earlier studies also showed that analytic models for arc cross 
sections are notoriously problematic and inadequate for an accu- 
rate absolute quantification of arc probabilities. However, since 
it was our main interest in this paper to quantify changes in 
the arc optical depth relative to cosmological-constant models, 
we applied reasonably flexible analytic lens models which we 
demonstrate to reproduce the relevant features of fully numeri - 
cal results (e.g. lBartelmann et alll998tlMeneghetti et all2003l) . 
as described below. 



We have chosen to use lenses with NFW density profile whose 
lensing potential was elliptically distorted. This model has sev- 
eral advantages. First, it naturally incorporates the dependence 
of central halo concentration on cosmology, through the halo 
formation time. Second, it agrees with density profiles consis- 
tently found in numerical simulations. Third, the elliptical dis- 
tortion strongly decreases th e mismatch between a nalytic and 
fully numerical lens models (Meneghet ti et all l2003). Although 
elliptical lensing potentials can lead to locally negative surface 
mass densities for sufficiently large ellipticities, this is a mild 
effect which happens only well outside the core for our model, 
thus it should be irrelevant for our purposes. We further adopted 
a simple estimate for the arc cross section and verified that this 
estimate scaled with cluster mass and dark energy parameter w 
in the same way as fully numerically determined arc cross sec- 
tions do. In particular, our simple analytic description for the arc 
cross section well reproduces the approximately quadratic scal- 
ing with cluster mass obtained fully numerically from the same 
elliptically distorted NFW cluster mass models. Although we 
would not trust this model for any accurate quantitative predic- 
tion of arc probabilities, we are confident that it can be used for 
the relative statements intended here. 

Our main findings are as follows: 

• Our simple lens model is indeed capable of reproducing the 
order-of-magnitude difference between ACDM and open 
CDM models found in the earlier numerical study. This 
shows that the change in halo concentration between the 
two models can explain the sensitivity of arc statistics for 
the cosmological constant. An additional effect into the 
same direction is contributed by the steep mass function of 
galaxy clusters. An increase in halo concentration lowers 
the minimum mass required for significant strong lensing, 
and this makes many more halos available for arc forma- 
tion. 

• Although increasing the dark energy parameter w has a no- 
ticeable effect on the optical depth for arc formation, it can- 
not increase the arc optical depth to a level compatible with 
that found in open CDM models. This result arises due to 
a combination of three main effects: first, individual halos 
of fixed mass get more concentrated in dark energy than 
in cosmological constant models; second, lower mass ha- 
los than before become able to form large arcs; and third, 
the requirement that arcs be large imposes a lower limit 
on cluster masses. Our analytic model thus suggests that 
arc statistics cannot be reconciled with low-density, spa- 
tially flat cosmological models which are now dominated 
by dark energy, i.e. the discrepancy between arc statistics 
and the cosmological model favoured by most, if not all, 
recent cosmological experiments is not expected to disap- 
pear if the dark energy is not a cosmological constant, but 
has an equation of state p = wpc 2 with w > — 1. 

Being based on several simple analytic estimates, this study 
can only provide a tentative answer. Detailed numerical simu- 
lations will be necessary for reliable absolute quantifications of 
the arc optical depth expected in cosmological models with dark 
energy instead of a cosmological constant; similarly, it is neces- 
sary to quantify the dependence of the effect we find here on the 
specific dark energy model considered, such as a quintessence 
scalar field. 

According to our preliminary analytic results presented here, 
it appears that the solution to the a rc statistics problem as de- 
scribed by Bartelmann et al. ( 1998) is probably not to be found 
in the cosmological model alone, but more in the details of clus- 
ter structure and the history of cluster evolution. The recent sug- 



gestion by Wambs ganss et alJ (1200 3) that it could simply be re- 
moved by considering higher-redshift sources is certainly not 
a viable solution because the problem arose from comparing 
simulations to a well-defined sample of highly X-ray luminous 
clusters which wer e observed for arcs down to a well-define d 
photometric limit (Le Fevre et all 1 19941 [Luppino et al., 1999). 
In combination with other recent resu lts on the strong- l ensing 
properties of dark-matte r halos (e.g. Ogu ri et all 1200 It [baud 
120021 lOeuri et al J 12003) and the puzzling discoveri es of very 
high-redshift arcs in high-redshift clusters ( Zaritskv & Gonzalez 
2003; Gladder s' et al.ll2003l) . this seems to imply that arc statis- 
tics are teaching us that the properties of individual clusters and 
their evolution over time is still insufficiently understood. 



Acknowledgements 



We wish to thank Peter Schneider, Achim Weiss and Simon 
White for useful discussions and comments. This work has 
been partially supported by Italian MIUR (Grant 2001, prot. 
2001028932, "Clusters and groups of galaxies: the interplay of 
dark and baryonic matter"), CNR and ASI. 



REFERENCES 



Abroe, M., Balbi, A., Borrill, J., Bunn, E., et al. 2002, MNRAS, 334, 11 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 
Bartelmann, M. 1996, A&A, 313, 697 

Bartelmann, M., Ehlers, J., & Schneider, P. 1993, A&A, 280, 351 
Bartelmann, M., Huss, A., Colberg, J., Jenkins, A., & Pearce, F. 1998, A&A, 
330, 1 

Bartelmann, M., Perrotta, F., & Baccigalupi, C. 2002, A&A, 396, 21 

Bartelmann, M., Steinmetz, M., & Weiss, A. 1995, A&A, 297, 1 

Bartelmann, M. & Weiss, A. 1994, A&A, 287, 1 

Bullock, J., Kolatt, T., Sigad, Y., et al. 2001, MNRAS, 321, 559 

Carroll, S. 2001, Living Reviews, 4, 1 

Chiba. M. & Yoshii, Y. 1999, ApJ, 510, 42 

Cooray, A. 1999, ApJ, 524, 504 

Douspis, M., Riazuelo, A., Zolnierowski, Y, & Blanchard, A. 2003, A&A. 405, 
409 

Efstathiou, G., Moody, S., Peacock, J., Percival, W., et al. 2002, MNRAS, 330, 
29 

Eke, V., Navarro, J., & Steinmetz, M. 2001, ApJ, 554, 1 14 
Falco, E., Kochanek, C, & Mufioz, J. 1998, ApJ, 494, 47 
Flores, R., Mailer, A., & Primack, J. 2000, ApJ, 535, 555 
Flores, R. & Primack, J. 1994, ApJ, 427, LI 
Gioia, I. & Luppino, G. 1994, ApJS, 94, 583 

Gladders, M., Hoekstra, H., Yee, H., Hall, P., & Barrientos, L. 2003, ApJ sub- 
mitted; preprint astro-ph/0303341 
Golse, G. & Kneib, J.-P. 2002, A&A, 390, 821 

Halverson, N, Leitch, E., Pryke, C, Kovac, J., et al. 2002, ApJ, 568, 38 

Helbig, P. 1999, A&A, 350, 1 

Kassiola, A. & Kovner, I. 1993, ApJ, 417, 450 

Kaufmann, R. & Straumann, N. 2000, Annalen der Physik, 9, 384 

Keeton, C. 2002, ApJ, 575, LI 

Kochanek, C. 1996, ApJ, 466, 638 

Kochanek, C. & White, M. 2001, ApJ, 559, 531 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

— . 1994, MNRAS, 271, 676 

Le Fevre, O., Hammer, F., Angonin, M., Gioia, I., & Luppino, G. 1994. ApJ. 
422, L5 

Lee, A., Ade, P., Balbi, A., Bock, J., et al. 2001, ApJ, 561, LI 
Lehar, J., Falco, E., Kochanek, C, et al. 2000, 536, 584 

Luppino, G, Gioia, I., Hammer, R, Le Fevre, O., & Annis, J. 1999, A&AS, 136, 
117 

Mainini, R., Maccio, A., & Bonometto, S. 2003, New A, 8, 173 
Melchiorri, A., Mersini, L., Odman, C, & Trodden, M. 2003, preprint astro- 
ph/0211522 

Meneghetti, M., Bartelmann, M., & Moscardini, L. 2003, MNRAS, 340, 105 
Meneghetti, M., Bolzonella, M., Bartelmann, M., Moscardini, L., & Tormen, G. 

2000, MNRAS, 314, 338 
Meneghetti, M., Yoshida, N., Bartelmann, M., et al. 2001, MNRAS,, 325, 435 
Navarro, J., Frenk, C, & White, S. 1996, ApJ, 462, 563 
— . 1997, ApJ, 490, 493 

Netterfield, C, Ade, P., Bock, J., Bond, J., et al. 2002, ApJ, 571, 604 
Oguri, M. 2002, ApJ, 573, 51 

Oguri, M., Lee, J., & Suto, Y. 2003, ApJ submitted; preprint astro-ph/0306102 
Oguri, M., Taruya, A., & Suto, Y. 2001, ApJ, 559, 572 
Peebles, P. & Ratra, B. 2002, Rev. Mod. Phys., 75, 599 

Perlmutter, S., Aldering, G, Goldhaber, G, Knop, R., et al. 1999, ApJ, 517, 565 
Porciani, C. & Madau, P. 2000, 532, 679 
Press, W. & Schechter, P. 1974, ApJ, 187, 425 



Quast, R. & Helbig, P. 1999, A&A, 344, 721 
Reiprich, T. & Bohringer, H. 2002, ApJ, 567, 716 
Richstone, D., Loeb, A., & Turner, E. L. 1992, ApJ, 393, 477 
Riess, A. G, Filippenko, A. V, Challis, P., et al. 1998, AJ, 116, 1009 
Sheth, R. & Tormen, G. 1999, MNRAS, 308, 119 

Smoot, G. 1999, in AIP Conference Proceedings, Vol. 476, 3K Cosmology, ed. 

L. Maiani, F. Melchiorri, & N. Vittorio 
Spergel, D., Verde, L., Peiris, H., Komatsu, E., et al. 2003, ApJ submitted; 

preprint astro-ph/0302209 
Wambsganss, J., Bode, P., & Ostriker, J. 2003, ApJ submitted; preprint astro- 

ph/0306088 
Wang, L. & Steinhardt, P. 1998, ApJ, 508, 483 

Wang, X., Tegmark, M., & Zaldarriaga, M. 2002, Phys. Rev. D, 65, 123001 
Wu, X.-P. & Mao, S. 1996, ApJ, 463, 404 
Zaritsky, D. & Gonzalez, A. 2003, ApJ, 584, 691 



